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1 Introduction 

Comparative statics is a special case of dynamics, in which a unique stable 
equilibrium is assumed to exist, and two equilibria — determined by different 
parameter values — are compared. Similarly, dynamics is a special case of 
an evolutionary process, in which the degrees of freedom of the system are held 
constant. Marshall, whose _Prmcip/es[^6| played a large role in setting economics 
upon the comparative statics course, was nonetheless aware of the primacy of 
the evolutionary paradigm when he wrote on page xiv that: 

The Mecca of the economist lies in economic biology rather than in 
economic dynamics. But biological conceptions are more complex 
than those of dynamics; a volume on Foundations must therefore 
give a relatively large place to mechanical analogies; and frequent 
use of the term 'equilibrium', which suggest something of a statical 
analogy... .The modern mathematician is familiar with the notion 
that dynamics includes statics. If he can solve a problem dynam- 
ically, he seldom cares to solve it statically also. . . But the statical 
solution has claims of its own. It is simpler than the dynamical; 
it may afford useful preparation and training for the more difhcult 
dynamical solution; and it may be the first step towards a provi- 
sional and partial solution in problems so complex that a complete 
dynamical solution is beyond our attainment. 

One century later, intellectual progress has enabled us to contemplate the 
construction of evolutionary models of processes which could once only be mod- 
eled statically. 

Early ecological theories followed a similar path of development. Interactions 
between organisms were studied in order to determine the population densities 
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that kept the ecology in equihbrium and the only form of dynamics considered 
was instantaneous movement from one equilibrium to another. It was always 
assumed that an ecology was always near a stable equilibrium, in spite of counter 
examples such as the stable limit cycles in the Lotka-Volterra predator-prey 
system Q. 

In the 1960s, interest in dynamical systems was reawakened, and the study of 
Chaos was born. These techniques were applied to ecology in the 1970s, and it 
was realised that the situation was vastly more complicated, with not only limit 
cycles being possible, but also entities called strange attractors. Over the years, 
appreciation has grown that the stable equilibrium is more of an exception than 
the norm. 

In the mid 1980s, Raup and Sepkowskip^ noticed certain statistical regu- 
larities with the pattern of speciation and extinction in the fossil record. These 
were later shown to be in the form of a power law[?, ?], i.e the pattern of extinc- 
tions was such that the frequency of an extinction event of size x is proportional 
to x~", where a is a small positive constant. Power laws crop up in many 
different areas, psf 

Per Bak[^ studied a model of a sandpile that had a continuous stream 
of sand added to the top. Bak's process goes by the name of self-organised 
criticality, and he, along with KauffmanpSf promote this as an explanation for 
Raup's data. It now seems likely that evolution is an endogenously self-organised 
critical process[] 

In my work, which I will introduce later, I have demonstrated the existence 
of self-organised criticality in a model evolutionary ecology. The criticality is 
quite a robust feature over a wide range of input parameters. This indicates 
that that evolutionary systems are typically endogenously critical. 

In summary, we have a hierarchy of approaches, from the static, to the 
dynamic to the evolutionary: 

statics C dynamics C evolution 

The state of ecological thought has followed this chain from left to right, as 
computational techniques have improved to embrace the computationally more 
difficult dynamical models, and then the even more computationally difficult 
evolutionary models. It is to be hoped that the same path will be followed by 
economic modeling. 

There is a clear parallel between the development of ecological thought and 
that of economic thought. Dynamics was introduced to economics and cham- 
pioned by Kaldor, Goodwin et al. in the 1950s and 60s; and boosted by the 
developments in nonlinear analysis in the 1980s, through the work of people 

^Alternatively, perhaps exogenous influences (such as volcanos, meteorites, climatic fluctu- 
ations, supernovae etc) obey a power law spectrum. It is quite possible that both endogenous 
and exogenous effects contribute to the statistical properties of biological evolution and that 
a future research program will look at untangling the two effects. To this effect, Newman[?] 
provides an interesting perspective. 
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like Blatt|^. This volume testifies to a growing industry of dynamical economic 
modeling. However, perhaps now is the time to embark on the next rung up 
the ladder, and approach economics from an evolutionary point of view. Then 
perhaps we might be catching a glimpse of the Mecca that Marshall referred to 
a hundred years ago. 

This paper first considers a general dynamical system undergoing evolution 
of the determining equation, then outlines Ecolab, a model of species interaction 
undergoing evolution, finally introducing a possible economics model based on 
the von Neumann model with evolution of the interaction coefficients. 



2 Linearisation of a Dynamical System 

Our launchpad for a theory of evolutionary systems is dynamical systems theory || 
Typically this will be manifested in a first order nonlinear differential equation 
of the form 

X = f (x) 

where x G M" and f : M — > M. The dot refers to the derivative with respect to 
time. 

Dynamical systems theory starts by considering the equilibria of the system, 
i.e. the points x such that f(x) — 0. Then in the neighbourhood of x, the 
behaviour of the system is determined by the linear approximation 

i-i?f|i-(x-ic). 

The stability of x is determined by the negative definiteness of I^f 1x0- This 
condition imposes n inequalities on the system constraining the form of f . There 
may additionally be a further n inequalities for x to be a meaningful solution, 
for example if the components of x are production values, every component of 
x must be non-negative. 

General equilibrium economics has attempted to find the conditions under 
which a unique, stable equilibrium will exist. This needn't be the case, and in- 
teresting (i.e. bounded) behaviour can take place around unstable equilibria, in 
the form of limit cycles or even the strange attractors beloved of chaos theorists. 
A favourite model of the latter researchers is the logistic equation, which first 
arose in a biological context but has been applied to economics [sT] amongst 
other things. 

That limit cycles and chaotic behaviour can be observed in economics is a 
view that should have by now been accepted[^. However, the question remains 
as to whether this behaviour is pathological, i.e. whether linear neoclassical the- 
ory is applicable in most cases, and the only remaining difficulty is determining 
if linear theory applies to a specific economy (the "econometric problem" ) , or 
whether chaotic behaviour is indeed the norm. 

■^Uf lii is negative definite if x ■ Df |x ■ x < for all x. This also implies that all eigenvalues 
of D{ have negative real part. 
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3 Limits of Linear Economic Theory 



Returning to the question of stability, the fact that 2n inequahties must be satis- 
fied would imply that a randomly chosen rt-dimensional economics model would 
have a probability of 4^" of a given equilibrium being stable. The situation 
does look bleak, but economics are not generated randomly in the real world, 
rather they are the result of an evolutionary process. We need to examine the 
process of cultural evolution to answer this question. 

The analogue of mutation in biology would be innovation in economics, as 
a new process or technique introduced into production, or as a new form of 
marketing, or a new company with a somewhat unusual approach to doing 
business. The effect of innovation is to add new degrees of freedom to the 
dynamical system, which usually will destabilise the system. The system will 
then tend to evolve so as to lose some of the degrees of freedom, by for example 
old production techniques being abandoned, or companies going bankrupt. This 
is analogous to species becoming extinct in the natural world. 

Lets consider what happens to the largest eigenvalue of Di\^. Suppose 
initially, the system has a stable equilibrium, in which case all the eigenvalues 
have negative real part. As innovations are added to the system, the largest 
eigenvalue will increase towards zero. As it passes zero, the system destabilises, 
and the system will start to exhibit limit cycles or chaotic behaviour. As further 
innovations are added to the system, a property called permanenc^ is no longer 
satisfied, and some event such as a bankruptcy will occur to remove active 
processes from the system. This will restore permanency to the system, and 
possibly even stability. Such a process is called self-organised criticality^ which 
gives rise to a power law spectrum of the booms and busts, successful innovations 
and bankruptcies.^ 



4 Ecolab and the Dynamics of Evolution 

This section outlines a model of an evolving ecology[^ that is analogous 
to an economic system with input-output relations of production and product 
innovation. The ecology is described by a generalised Lotka-Volterra equation, 
which is perhaps the simplest ecological model to use. 

■hi = TiTli + ^ PijUiUj (1) 
J = l 

•^that there is a set of points xq whose trajectories x(t) always remain away from the 
boundary, i.e. Xi{t) > S 3<5 > 

*This really implies that conclusion derived from comparative static economic analysis 
are almost never valid, except perhaps on sufficiently small time scales while the maximum 
eigenvalue of Df |x is negative. 
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Here r is the difference between the birth rate and death rate for each species, 
in the absence of competition or symbiosis. (3 is the interaction term between 
species, with the diagonal terms referring to the species' self limitation, which 
is related in a simple way to the carrying capacity Ki for that species in the 
environment by Ki = —riPa. In the literature (eg StrobeckQ, CaseQ) the 
interaction terms are expressed in a normalised form, a.^ = ~Ki/ri(3ij, and 
an = 1 by definition, n is the species density. 

These equations are simulated on a simulator called Ecolab. ||2^ The vectors n 
and r are stored as dynamic arrays, the size of which (i.e. the system dimension) 
can change in time. The interaction array is stored in row/column sparse form, 
consisting of the four arrays diag, val, row and col. Equation (|l|) can be 
written as: 

tmp[row] = beta. val * n[beta.col]; 
n += (r + beta. diag + tmp) * n; 

This code makes up the generate operator in the Ecolab system. Other op- 
erators include compact, which removes species that have become extinct from 
the system (to optimise computational performance) and mutate, which adds a 
certain number of new species to the system, according to a specific algorithm to 
be discussed later. The operators can be called from a scripting language called 
TCL|21|, that allows different types of experiments to be performed without 
recompiling the code. 

Before discussing the mutation algorithm in more detail, equation (|l|) must 
be analysed to determine the conditions /3 must satisfy for the system to be 
real, and also to determine the different regimes of dynamics, from the linear 
(stable equilibrium) case, to limit cycles and chaos to the actual breakdown of 
the ecosystem. 



4.1 Linear Analysis 

Linear analysis starts with the fixed point of equation (|l|) 

n = -l3-'r, (2) 

where n = 0. There is precisely one fixed point in the interior of the space of 
population densities (i.e. n such that > 0) provided that all components of 
n are positive, giving rise to the following inequalities: 

n,= {f3-'r)^>0, (3) 

This interior space is denoted R^^^ mathematically. 

There may also be fixed points on the boundary of where one or 

more components of n are zero (corresponding to an extinct species). This is 
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because the subecology with the Uving species only (i.e. with the extinct species 
removed) is equivalent to the full system. 

The stability of this point is related to the negative definiteness of derivative 
of h at n. The components of the derivative are given by 



Substituting eq (||) gives 



dfij 
djij 



(4) 



(5) 



Stability of the fixed point requires that this matrix should be negative 
definite. Since the (/3^^r)^ are all negative by virtue of this is equivalent 
to (3 being negative definite, or equivalently, that its rigp eigenvalues all have 
negative real part. Taken together with the inequalities this implies that 
2nsp inequalities must be satisfied for the fixed point to be stable. This point was 
made by Strobeck|^, in a slightly different form. (Note that Strobeck implicitly 
assumes that rifii/ Ki > 0, so comes to the conclusion that 2nsp— f conditions 
are required.) If one were to randomly pick coefficients for a Lotka-Volterra 
system, then it has a probability of 4^"'>'p of being stable, i.e. one expects 
ecosystems to become more unstable as the number of species increases . 



4.2 Permanence 

Whilst stability is a nice mathematical property, it has rather less relevance 

when it comes to real ecologies. For example the traditional predator-prey 

system studied by Lotka and Volterra has a limit cycle. The fixed point is 

decidedly unstable, yet the ecology is permanent in the sense that both species' 

densities are larger than some threshold value for all time. Hofbauer et al. |^ 

and Law and Blackfordpq] discuss the concept of permanence in Lotka-Volterra 

— n 
systems, which is the property that there is a compact absorbing set JV[ c M^''^ 

i.e once a trajectory of the system has entered M., it remains in M. They derive 

a sufficient condition for permanence due to Jansen[Tl}| of the form: 

^Pifiins) = ^P^{r^ - ^PrjfiBj) > 0, 3p, > (6) 

i i j 

for every fiB equilibrium points lying on the boundary (nsi = 3«), provided 
the system is bounded (or equivalently dissipative). This condition is more gen- 
eral than stability of the equilibrium — the latter condition implies that a 
local neighbourhood of the equilibrium is an absorbing set. Also, the averaging 
property of Lotka-Volterra systems implies that the equilibrium must lie in the 
positive cone So (^) must still hold for permanence. 
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Consider the boundary points fiB that are missing a single species i. Then 
Jansen's condition for these boundary points is 

"Ti - '^PijfiBj > 0. (7) 

j 

This set of conditions is hnearly independent. Let the number of such boundary 
points be denoted by < risp. Then the set of conditions wiU have rank 
ns V < Ugp (the number of hnearly independent conditions) , so the system 
has at most probability 2~^^p~'^ of satisfy Jansen's permanence condition if the 
coefficients are chosen uniformly at random. As stability is also sufficient for 
permanence, the probability lies between 4^""*^? and 2""'^?^'^. 

Another rather important property is resistance to invasion.^ Consider a 
boundary equilibrium fiB- If it is proof against invasion from the missing species, 
then the full system cannot be permanent. For the boundary points that miss a 
single species, this implies that condition is necessarily satisfied for perma- 
nence, along with (^). The probability of permanence is then bounded above 
by 2""'>5P""«. 

Thus whilst a randomly selected ecology is more likely to be permanent than 
to have a stable equilibrium, the likelihood decreases exponentially with increase 
in species number. 



4.3 Boundedness 

It is necessary that the ecology be bounded, ie that < 3N, Vt > 0. 

This requires 

^li, ^ r ■ n + n ■ f3n<0, yn -.^72, > N 3N (8) 

i i 

As n becomes large in any direction, this functional is dominated by the 
quadratic term, so this implies that 

n • /3n < Vn : > 0. (9) 

If strict equality holds, then r ■ n < 0. Negative definiteness of (3 is sufficient, 
but not necessary for this condition. Another sufficient condition is to require 
< and Pij + (3ji < 0, which is used in the current study. This 
condition is satisfied by the Predator-Prey equations, and so does allow multi- 
trophic systems to be built, but does not allow the possibility of symbiosis. 
Its main advantage is its simplicity of implementation, along with the range of 
interesting (i.e. non limit point) behaviour it encompasses. 

4.4 Mutation 

Adding mutation involves adding an additional operator to equation (^ 

ri = r*n + n* (3n + mutate(/i,, r, n) (10) 
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where * refers to elementwise multiplication. 

The mutation operator must generate new degrees of freedom i > Tigp (where 
Usp is the number species currently in the ecology), somehow defining the new 
ecological coefficients {ri\i > nsp},{l3ij\i > Usp or j > Ugp} from the previous 
state of the system. In reality, there is another layer (hidden in equation (^ 
called the genotypic layer, where each organism has a definite genotype. There 
is a specific map from the genotypic layer to the space of ecological coefhcients 
(hereafter called the phenotypic layer) called the embryology. Then the mutation 
operator is a convolution of the genetic algorithm operations operating at the 
genotypic layer, with the embryology. 

A few studies, including Ray's Tierra world, do this with an explicit mapping 
from the genotype to to some particular organism property (e.g. interpreted as 
machine language instructions, or as weight in a neural net). These organisms 
then interact with one another to determine the population dynamics. In this 
model, however, we are doing away with the organismal layer, and so an ex- 
plicit embryology is impossible. The only possibility left is to use a statistical 
model of embryology. The mapping between genotype space and the popula- 
tion parameters r, (3 is expected to look like a rugged landscape, however, if two 
genotypes are close together (in a Hamming sense) then one might expect that 
the phenotypes are likely to be similar, as would the population parameters. 
This I call random embryology with locality. Here, we tend to idealise genotypes 
as bit strings, although strings over an arbitrary alphabet (eg the four DNA 
bases ACGT) can equally be considered. |^ 

In the simple case of point mutations, the probability P{x) of any child lying 
distance x in genotype space from its parent follows a Poisson distribution, as 
this is the distribution of the number of bit flips, or deletions that might occur 
with a point mutation. Random embryology with locality implies that the phe- 
notypic parameters are distributed randomly about the parent species, with a 
standard deviation that depends monotonically on the genotypic displacement. 
The simplest such model is to distribute the phenotypic parameters in a Gaus- 
sian fashion about the parent's values, with standard deviation proportional 
to the genotypic displacement. This constant of proportionality can be con- 
flated with the species' intrinsic mutation rate, to give rise another phenotypic 
parameter fi. It is assumed that the probability of a mutation generating a pre- 
viously existing species is negligible, and can be ignored. We also need another 
arbitrary parameter p, "species radius" , which can be understood as the min- 
imum genotypic distance separating species, conflated with the same constant 
of proportionality as fi. 

We may represent the Ecolab embryology as a probability distribution /(p, g) = 

I2.B^_JI!_^ where p ^ \r, ~ rj\/\r,\ or p = \(3ik - Pjk\/\f3ik\ is the distance 




between two species' phenotypic parameters, and g is the difference between the 

^The Hamming distance is the number of bits (bases) that differ between the two strings. 
So for example if a single bit has been removed from one string, the Hamming distance is one. 
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Figure 1: Probability distribution of the relation between genotype difference 
and the corresponding phenotype difference 

two genotypes. Figure |l| shows the general form of this probability distribution. 

Figure |^ shows the probability distribution of a mutant phenotypical coeffi- 
cient about that of its parent's value. This is given by 



In summary, the mutation algorithm is as follows: 

1. The number of mutant species arising from species i within a timestep is 
^^ifi'r^i/ p. This number is rounded stochastically to the nearest integer, 
e.g. 0.25 is rounded up to 1 25% of the time and down to 75% of the 
time. 

2. Roll a random number from a Poisson distribution e~^li^^P to determine 
the standard deviation a of phenotypic variation. 

3. Vary r according to a Gaussian distribution about the parents' values, 
with otq as the standard deviation, where tq is the range of values that r 
is initialised to, ie rg = max^rilj^g — minirijj^Q. 




(11) 
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p 

Figure 2: The probability distribution of a mutant pheno typical coefficient 
about that of its parent's value. This has been calculated by numerical in- 
tegration from formula (lO). Note that the curve actually diverges at 0. 
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4. The diagonal part of (3 must be negative, so vary /3 according to a log- 
normal distribution. This means that if the old value is /3, the new value 
becomes /?' = — exp(— ln(/3) + a). These values cannot arbitrarily ap- 
proach 0, however, as this would imply that some species make arbitrarily 
small demands on the environment, and will become infinite in number. 
In ecolab, the diagonal interactions terms prevented from becoming larger 
than — r/(.l * INT_MAX), where r is the corresponding growth rate for the 
new species. 

5. The off diagonal components of /3, are varied in a similar fashion to r. 
However new connections are added, or old ones removed according to 
[l/pj, where p £ (—2,2) is chosen from a uniform distribution. The 
values on the new connections are chosen from the same initial distribu- 
tion that the off diagonal values where originally set with, ie the range 
miui^j Ij^Q to maxi^j /3y Ij^p. Since condition is computationally 
expensive, we use a slightly stronger criterion that is sufficient, compu- 
tationally tractable yet still allows "interesting" non-definite matrix be- 
haviour namely that the sum -I- (3ji should be non positive. 

6. fi must be positive, so should evolve according to the log-normal distribu- 
tion like the diagonal components of /3. Similar to (3, it is a catastrophe 
to allow fj, to become arbitrarily large. In the real world, mutation nor- 
mally exists at some fixed background rate — species can reduce the level 
of mutation by improving their genetic repair algorithms. In ecolab, this 
ceiling on fi is given by the mutation(random,maxval) variable. 

4.5 Typical Results 

Figure]^ shows the time behaviour for the number of species in the ecosystem for 
a typical run. The phenotypic parameters were seeded randomly in the ranges 
-0.005 < r < 0.01, -5 x 10"^ < /3diag < -1 x 10-^ -0.001 < /S^ffdiag < 
0.001 and < /x < 0.09. The r and f3 values were chosen so that several 
hundred individuals will be supported in the case of a single species system, and 
the offdiagonal terms large enough to permit interesting interactions between 
species, but not so large that the system collapsed to zero immediately, p was 
set at 10^, which was chosen by examining the histogram of differences between 
all the species. If p was too small, then a species' mutant offspring would be 
too similar to its parent to be really a new species. This shows up as a peak at 
small separation values of the histogram, which shouldn't be there according to 
the law of competitive exclusion. 

The system rapidly evolves to one of the fixed points (by a massive extinction 
event!) with a negative definite /3. Over time, mutations build up in the system, 
decreasing the stability of the system. What then follows are periods of episodic 
extinctions, and system growth through speciation. This is an example of self 
organised criticality^, and gives rise to power law behaviour. 
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Figure 4: Distribution of species lifetimes on a log-linear plot. Distribution is 
unnormalised. Horizontal scale is the natural logarithm of species lifetime in 
timesteps. 

Do we see the same power law behaviour observed by others?^ |l| The an- 
swer is emphatically yes. If speciation and extinction events occurred uniformly 
throughout history, as one might naively expect, one would expect a Poisson 
distribution for species lifetimes. On a log-linear plot, this would be a straight 
line. Alternatively, if a power law spectrum was evident, the log- log plot would 
be straight. The two plots are shown in figures ^ and ^. Effectively, this is 
telling us that not only is there not a stable ecological equilibrium, there isn't 
even a steady state, whereby extinctions are balanced by speciation (a common 
ecological assumption) . 

This model then is a concrete example of the self-organised criticality pre- 
dicted in these types of systems in section??. The next section examines a 
possible economic model that is analogous to Ecolab, and could even be im- 
plemented using the same simulation software. It would be surprising if the 
dynamics weren't critically self-organised. 
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5 Building an Economic Dynamics 



Many inferential similarities can be drawn between the biological evolutionary 
model of Ecolab and the processes of a capitalist economy. The obvious analogy 
for a biological species is a product, and for Darwinian evolution the process 
of technological change. I consider a model economics (Econolab) based on 
the insights of von Neumann, one of the founders of complexity theory, who 
introduced von Neumann Technology in the late 1930s |]3^ [33| . In this model 
economy, there is a set of commodities labeled i E Af = {1 . . . N}), and a set 
of technologies or processes labeled m G M — {1...A/}. Each process has 
an activity z™, input coefficients ami and output coefficients bmi, such that in 
one time step, amiZm of commodity i (amongst others) is consumed to produce 
bjnjZj of commodity j (amongst others). The coefficients ami and b^i may be 
zero for some values of m and i, corresponding respectively to processes that 
do not require a particular input, or do not produce a particular output. This 
differs from von Neumann's original approach, and is more in line with that of 
Kemeny, Morgenstern and Thompson Jl^. Blatt[^ gives a good introduction to 
this model, discussing it flexibility in dealing with a range of economic processes. 
In the words of Blatt (p67) : 

The von Neumann work is a great achievement of mathematical 
model building in dynamic economics. It is the best available theory 
of capital and of rate of return. 

That said, there are many issues of significance in capitalism which are not 
captured in the von Neumann method, and which cannot be modeled in an initial 
rendition of Econolab. These include effective demand[?], income distribution, 
variable capacity and utilisation, credit and debt[?, ?]. 

To relate von Neumann's work back to the Ecolab ecological model, the 
input/output coefficients ami/bmi are fixed like the r^, pij or equation (|l]), and 
Zm is a free variable like n^. In von Neumann's work, the dynamics is imposed 
in the form of an exponential growth condition: 

z„(t + l) =az,„(t) Vme (12) 

However, rather than assuming a particular form for the dynamics, we should 
be looking for a first order differential equation (or its difference equation equiv- 
alent) that describes the dynamics. Consider the monetary value of capital Km 
associated with process m. The rate of change of this capital may be written: 

(N N \ 

^ ^ ^7niPi ^ ^ ^miPi j ; (-^*^) 
1=1 1=1 / 

where pi is the price of commodity i. This has introduced two new sets of 
free variables Km and pi, for which we need to find closure relations. Clearly, 



15 



activity is limited by the availability of capital (we do not allow the possibility 
of credit here): 



N 



'^amiViZm < K,n (14) 

1=1 

For simplicity, let us assume that each process invests a fixed proportion of 
its capital into production, i.e. 

N 

^ ^ ImiPiZm — '^m-^mj '■ < Kjn ^ 1 (1^) 

i=l 



Substituting (H) into (|T|) 



"-m I j^j -L 



(16) 



This then, is a model dynamics analogous to the Lotka-Volterra equation (|i|). 
If price is a fixed quantity (as assumed in von Neumann theory) then ( |l6| ) is 
equivalent to the ansatz (|l2[). This is the equilibrium situation, rather like 
assuming that n — h. 

In reality, prices arc not fixed, and must have their own dynamics. The 
simplest way to do this is to look for a closure relation, that relates prices to 
activities. The neoclassical and Austrian traditions propose that price dynamics 
should act as a negative feedback on the activity dynamics (eq. (p^)), whereas 
the P-K and Sraffian tradition do not see price as an equilibrating mechanism 
[?]. In this work, however, we propose an ansatz on the form of the negative 
feedback, in a similar fashion to the ansatz used by Nose and Hoover to 
describe the thermostat that regulates the temperature of a non-equilibrium 
steady state system in a heat bath: 

/demand \ / Em=i a™^™ i\ /1'7^ 

(i^-r-^'"(E;L j"" ' ' 

This differs from von Neumann, who assumes that demand never exceeds supply, 
and if supply exceeds demand (i.e. a surplus), then the commodity is free 
{pi = 0). This would imply pi = 0, freezing prices. In effect this makes the 
system very stiff — equation (^^ softens the dynamics with tTj controlling the 
stiffness. 



6 Adding Evolution 

Now that we have an economic dynamics established, we need to consider how 
to develop an analogy between ecological and economic evolution. By direct 
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analogy with Ecolab, it is clear that when a process exhausts its capital {Km = 
0), it forever remains that way, so this is equivalent to extinction in ecosystems. 
Adding new processes and commodities is conceptually easy. Blatt p57-58Q: 

What about technological progress? This can be included by assum- 
ing that the list of activities m = 1,2, ...,M is not final, but new 
activities may be invented and hence become available for use, as 
time goes on. This makes the total number of processes a function 
of time: M — M{t). . . Von Neumann himself developed his theory on 
the basis of an unchanged technology (all input cocfhcients, output 
coefficients and the number of processes M are constant in time), 
and his successors have done the same. The inclusion of techno- 
logical progress appears to us to be a highly interesting avenue for 
further exploration. 

The difficulty is deciding how to choose new coefficients a^i , i>mi i i^i and tt^ 
when a new process is added. There is no genotype of a process — the closest 
thing to it is Dawkins's meme, and there is no genetic algorithm theory of 
the meme. Clearly new processes arise evolutionarily, with the new processes 
modeled on the old. The new coefficients will be varied randomly about the old 
values according to some kind of central distribution. 

Recent results from Ecolab indicate that the emergent dynamics of the sys- 
tem is rather insensitive to the specific type of mutation algorithm chosen. Work 
is currently under way to classify exactly what effects different assumptions 
make. 

In 1962, Arrow Q pointed out that the cost per unit for production of an 
artifact falls as an inverse power of the number of units produced: 

cost/unit oc iV^" 

This power law is most likely a consequence of the dynamics of technologi- 
cal innovation, relating to the statistical properties of the underlying "fitness" 
landscape, as it can be seen in Kauffman NK model[^. Presumably an evo- 
lutionary algorithm that searches process (and commodity) space according to 
the same power law would be optimally matched to generating change, how- 
ever another search algorithm would probably generate the same distribution 
of successful innovations, albeit on a different temporal scale. It should also 
be pointed out that large changes of process are likely to cost proportionally 
more than smaller changes. As any research budget is finite, the distribution of 
process improvements must therefore be finitely integrable (have a finite area 
underneath the curve), which the power law distribution is not, but the normal 
(Gaussian) distribution is. 
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7 Conclusion 



Economics is clearly a dynamic process, which given its complexity will be 
poorly described by a linear approximation about a stable equilibrium. Rather 
the properties of the equilibrium will be determined by cultural evolution which 
operates over a longer timescale than economics. It is likely that cultural evo- 
lution will produce a self-organised critical system, and this would be one of 
the first questions to study. Other questions that might be looked at include 
looking for evidence for the Arrow law, and looking for analogues to various 
biological laws, such as the species- area law 0, and dependence of biodiversity 
with latitude. 

Perhaps the most important point I would like to make is that rather than 
studying a finite dimensional dynamical system, we should be studying what 
might be called "open dimensional" dynamical systems, where the number of 
degrees of freedom is finite, but not fixed at any point in time. These systems 
must lie between finite dimensional spaces and infinite dimensional "functional 
analysis" type spaces. Only then might we achieve Marshall's economic biology, 
and have an understanding of why economic systems have evolved to be the 
way they are. 
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